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Abstract. The eigenvalue spectrum of the fractional quantum harmonic oscillator is calculated numerically 
solving the fractional Schrodinger equation based on the Riemann and Caputo definition of a fractional 
derivative. The fractional approach allows a smooth transition between vibrational and rotational type 
spectra, which is shown to be an appropriate tool to analyze IR spectra of diatomic molecules. 
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1 Introduction 

Besides Newton's principia Huygens treatise on the pen- 
dulum clock pQ may be considered as one of the most in- 
fluential works in physics. Since then, the description of 
periodical motion in terms of the classical harmonic oscil- 
lator is a standard example for an exactly solvable problem 
in classical mechanics. 

In non-relativistic quantum mechanics an arbitrary po- 
tential may be approximated using a harmonic expansion 
around the equilibrium position. Therefore it is the ideal 
tool to model vibrational degrees of freedom in many dif- 
ferent quantum systems. 

Consequently one of the first applications of the quan- 
tum harmonic oscillator was the analytic description of 
the vibrational energy contribution in quantum mechani- 
cal models of diatomic molecules [2]. 

To establish the standard model for a successful de- 
scription of excitation spectra of diatomic molecules in 
addition to vibrations, rotational degrees of freedom have 
to be considered, so that the total energy of such a system 
is given in the simplest approximation mainly as a sum of 
both contributions (including e.g. anharmonic corrections 
and a rot-vib interaction term): 



vib 



E Ti 



(i) 



It is the basic intention of this article, to propose an alter- 
native approach based on fractional calculus, which over- 
comes the traditional distinction between rotational and 
vibrational degrees of freedom by introducing a more uni- 
versal description, which treats rotations and vibrations 
simultaneously within the generalized framework of frac- 
tional oscillations. 
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The fractional calculus [3] -[5] provides a set of axioms 
and methods to extend the standard derivative definition 
in a reasonable way from integer order n to arbitrary order 

^ - (2) 

and therefore allows for an extended, smooth derivative 
definition, which may be used to determine generalized 
symmetries, which go beyond the U(n)- symmetry of the 
standard harmonic oscillator. 

In the following, we will first present the main results of 
a numerical solution of the fractional Schrodinger equation 
with the fractional harmonic oscillator potential. 

We will then propose a fractional analogue to the stan- 
dard rot-vib model used for a description of spectra of di- 
atomic molecules. As a first application, we will use this 
model to describe the IR-spectrum of hydrogen chloride 
H 35 C1. 



2 The fractional derivative 

We will investigate the spectrum of the fractional quan- 
tum harmonic oscillator for two different definitions of the 
fractional derivative, namely the Riemann [6] and Caputo 
[7] fractional derivative. Both types are strongly related. 

Starting with the definition of the fractional Riemann 
integral 

R I a f(x)= (3) 

r (*/?/)(*) = / X di (x - o a -V(0 x > o 
\( R i«f)( X ) = T ^f°dn(z-xr- i m) x <o 

the fractional Riemann derivative is defined as the product 
of a fractional integration followed by an ordinary differ- 
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entiation: 



8 

£\Ot rl — a 

R x ~ dx R 



(4) 



The Caputo definition of a fractional derivative follows an 
inverted sequence of operations Q. An ordinary differen- 
tiation is followed by a fractional integration 



f) a — T 1 



dx 



(5) 



Applied to a function set f(x) — x na using the Riemann 
fractional derivative definition Q we obtain: 



j:> (J : j. .L 



r(l + na 



r(l + (n-l)a) 



1— x>0 (6) 



while for the Caputo definition of the fractional derivative 
it follows for the same function set: 



3 a? x % 



n>0 







n = .x > 



(7) 



Thus for e.g. polynomials with x na both derivative defi- 
nitions only differ in the case n = 0. 



3 Fractional quantum mechanics 

The transition from classical mechanics to quantum me- 
chanics may be interpreted as a transition from indepen- 
dent coordinate space and momentum space to a Hilbert 
space, in which space and momentum operators are treated 
similarly. 

Consequently one postulate of quantum mechanics sta- 
tes, that derived results must be independent of the spe- 
cific choice of e.g. a space or momentum representation. 
This is the mathematical manifestation of wave-particle- 
duality: A description in terms of either position x or wave 
vector k covers the properties of the same quantum object 
equivalently. 

This implies, that a successful fractional extension of 
quantum mechanics has to treat coordinates and conju- 
gated momenta equivalently: 

d n d a 
{x n ,-^}^K\^— } x>0 (8) 
1 ' dx n s 1 ' dx a J w 

This simultaneous treatment is the major difference be- 
tween a classical and a quantum mechanical treatment. 
Furthermore this is the key information needed to quan- 
tize any classical quantity like a Hamilton function. 

There are several different approaches, like canonical 
quantization [8] , [9] , path-integrals [10] , Weyl phase-space 
quantization [11] or stochastic quantization [T2] . 

To emphasize the equivalence between a space and a 
momentum representation, we will use the classical proce- 
dure, proposed by Dirac: The classical canonically conju- 
gated observables x and p are replaced by a pair of quan- 
tum mechanical observables which are introduced 



as derivative operators on a Hilbert space of square in- 
tegrable wave functions /. The space representations of 
these operators are: 



xf(x) = xf(x) 
pf(x) = -ihd x f(x) 



(9) 
(10) 



In order to extend this approach to the fractional case, we 
first make the following statement on parity conservation: 
The parity properties of the standard coordinate x and 
derivative d x should be conserved in the fractional case 
too. Therefore we postulate the following extensions valid 
on R: 



x a = sign(x)\x\ a 
D x = sign(» RjC <9£| 



(11) 
(12) 



and introduce the fractional pair {X,P} of conjugated 
fractional derivative operators, which are given in space 
representation: 



(l-a) 

Xf(x a )=[ — \ x a f(x a ) 



Pf(x a ) 



-) 

mc J 

h 



(13) 

mcD a x f{x a ) (14) 



The attached factors (h/mc)^ 1- ^ and (h/mc) a mc ensure 
correct length and momentum units. For the special case 
a = 1 these definitions correspond to the classical limits 
([9| and ([lOj). 

With these operators, we may quantize the classical 
Hamilton function of the harmonic oscillator: 



4 The fractional quantum harmonic oscillator 

The classical Hamilton function of the harmonic oscillator 
is given by 



tt P , 1 2 2 

tL class — — mcj X 

2m 2 



(15) 



Following the canonical quantization procedure we replace 
the classical observables {x,p} b y th e fractional derivative 
operators {X,P} according to (I13| and (14). The quan- 
tized Hamilton operator H a results: 



P 2 1 

H a = + -muj 2 X 2 

2m 2 

The stationary Schrodinger equation is given by 
H a & = 



(16) 



2m \mc J 



mc 



2(l-a) 

2a\ 



(17) 



Introducing the variable £ and the scaled energy E: 

r \—\ — ) ■>■•' (18) 



mu) j % 
V h \mc 
E' = huE 



(19) 
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Table 1. Parameters in the R;C ^ r ^ b (^,«) series expansion (24) 
and validity ranges of the fractional parameter a, which fulfills 
the requirement of normalizability of the wave function. 



type 



parity 



Riemann even 

Riemann odd 

Caputo even 

Caputo odd 



r 7i normalizable « 0) 

a - 1 0.5 < a < 2 o(£° 

a-1 a 0.25 <a< 2 o(f~ 

< a < 2 o(f°) 

a < a < 2 o(f a ) 



: 2c*-lx 



we obtain the stationary Schrodinger equation for the frac- 
tional harmonic oscillator in the canonical form 



1 



\e a W n (0 = E(n:a)V n (0 (20) 



In contrast to the standard quantum harmonic oscillator 
this fractional Schrodinger equation has not been solved 
analytically until now. 

An approximate solution for the energy levels has been 
derived by Laskin [13] within the framework of WKB- 
approximation, which is independent from a specific choice 
of a fractional derivative type: 



3 (n,a) = (rc+ -) 



T a/2 



l±a\ 
2a > 



n = 0,1,2,... 

(21) 

In a previous work [14] we have already emphasized, that 
these levels allow for a smooth transition from vibrational 
to rotational types of spectra, depending on the value of 
the fractional derivative coefficient a. 

£w KB (n, a * 1) ~ n + ^ n = 0, 1, 2, ... (22) 

£ WKB (n, a * 2) ~ (n + ^) 2 = n(n + 1) + 1/4 (23) 

Therefore the fractional harmonic oscillator may be a use- 
ful tool to describe rot-vib spectra, with a in the range of 
1 < a < 2. Consequently, the fractional harmonic oscilla- 
tor could be of similar importance as its classical counter- 
part, an exact solution is of fundamental interest. 

The main results from a numerical solution of the cano- 
nical Schrodinger equation of the fractional quantum har- 
monic oscillator (20) using the Riemann and Caputo frac- 
tional derivative definition are presented in the next sec- 
tion. 



5 Numerical solution 



For a numerical solution of (20), we expand the wave func- 



tion in a fractional Taylor-series with good parity 77: 

N 

a,c^(U)=ff^fl i (£)e 2ai £>0 (24) 



(25) 



where the parameters r and tt determining the type (Rie- 
mann or Caputo) of the fractional derivative definition 
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Fig. 1. Logarithmic plot of the first 10 energy levels of the 
fractional quantum harmonic oscillator in the range 0.85 < 
a < 2 based on the Riemann (blue), Caputo (red) and WKB- 
approximation (green) of the fractional derivative definition. 
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Fig. 2. Result of a least square fit of energy level spectra based 
on the Riemann (blue lines) and Caputo (red lines) derivative 
definition with the WKB-approximation ( |21| ). Plotted is the 
difference of a values for a given number of energy levels (n = 
10, 100) with fitted a WKB . 



used and parity of the wave function are listed in table 
[l] The ± index indicates the parity for n even and odd 
respectively. 

The coefficients (E) are determined using a standard 
shooting method: 

From the requirement, that the wave function should 
vanish at infinity for a given eigen value E n follows a de- 
termining condition for sufficiently large distance £ det : 



lim R>c 8£(f drt ,a,JS) = 



(26) 



which allows for an iterative procedure to calculate E n 
with arbitrary precision. 

For practical calculations, we determined the eigen val- 
ues with a precision of 32 significant digits with the set- 
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Fig. 3. Logarithmic probability density |^~ = n($, o>)\ 2 for the solution of the fractional quantum harmonic oscillator based on 
the Caputo (left) and Riemann (right) definition of the fractional derivative in the full range of allowed a values for real energy 
eigenvalues. The dashed line indicates the expectation value for the modified position operator (£). 



tings 20 < i de t < 35 and 1000 < N < 10000, depending 
on a. 

The major restriction for the range of allowed a values 
is the requirement of normalizability of the wavefunction 
obtained. From 



J oo 



dUc^t(^a) R , c ^^,a)* < oo 



(27) 



follows an upper bound a < 2. Of course, for a > 2 con- 
dition (26) may still be used as an equivalent of a box- 



normalization condition, but the results are no solutions 
of the harmonic oscillator potential any more. 

In figure [I] we compare the resulting energy level spec- 
tra to the WKB-approximation. While in the vicinity of 
a = 1 differences are negligible, for larger a values espe- 
cially for low lying energies the calculated spectra differ 
significantly, but as a first remarkable result we obtain 
a smooth transition from vibrational to rotational type 
of spectrum using the Riemann and Caputo fractional 
derivative definition respectively. 

The results of a fit of the energy levels with the WKB- 
approximation are plotted in figure [2] We may deduce, 
that for increasing number of energy levels n eq. (21) ap- 



proximates the energy values better, from this figure we 
may also deduce, that rotational spectra, which are char- 
acterized by 

^wkb — 2 are reproduced using the Caputo 
derivative for a slightly smaller a value, a c < 2, while for 
the Riemann derivative we have to use a slightly larger a 
value, a K > 2 and therefore using the Riemann fractional 
derivative definition implies the need for box normaliza- 
tion, if we want to generate the exact level spacing of 
rotational type. 

The behavior of the eigen functions differs 
significantly in the region 1 < a < 2 and a < 1. In fig- 
ure [3] we have plotted the occupation probability for the 
eigen function ^~ =11 (^a) for the Caputo and Riemann 
solution. 



Introducing the modified position operator £ 
its expectation value 



(0 



(28) 



yields the position information on the positive half-axis. 
This value is almost constant for a > 1. For increasing a 
only the side maxima are changing from inside to outside 
and therefore mark a smooth transition from Hermite-type 
to Bessel-type polynomials. 

In contrast for a < 1, the position value tends very 
fast to zero which means, that the wave function becomes 
strongly located at the origin in this case. 

As a consequence, the energy diagram for a < 1 is 
dominated by the kinetic term, since for the potential near 
the origin 

lim|£| 2a = (29) 

holds. 

This means, that the solutions of the fractional Schrodi- 
nger equation for the harmonic oscillator and the solutions 
for the free fractional Schrodinger equation more and more 
coincide for a < 1, a behavior, which we expect for any 
potential, which vanishes at the origin. 

The eigen functions of the free fractional Schrodinger 
equation are known analytically and are given in terms of 
the Mittag-Lemer functions E a (z) and E a ^(z) as: 



^+^(H,a)=£ a - 1 E 2a:a (-e a ) 



(30) 
(31) 
(32) 
(33) 



and the corresponding eigen values are determined in box- 
normalization from the zeros of these functions. For a < 1 
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Table 2. Parameter sets for the standard rot-vib-model (34) 
in the first row [2 and for the proposed generalized fractional 



model (139} to determine the energy levels of H 35 C1. Finally the 
rms-error is given. 



stan- 
dard 



10.59341 0.30718 

flUJe TlUJeXe 

2990.946 52.8186 



D e 

5.3194* 10" 
0.2243 



Pe 

7.51 * 10"' 

rms-error 

3.8% 



WKB 


B e 


a e 


9o 


TiUJe 




7.406 


0.202 


152.016 


2707.77 




OLJ 






rms-error 




1.997 


1.051 




3.3% 


Cap 


Be 




9o 


TlUJe 




5.264 


0.155 


-64.973 


2938.77 






Oiv 




rms-error 




1.91 


0.914 




3.8% 


Rie 


Be 




9o 


TlUJe 




2.588 


0.080 


10.355 


3095.36 




aj 






rms-error 




2.081 


0.937 




3.8% 



there exists only a finite number of zeros and therefore the 
energy spectrum is limited [T5] . 

In the upper row of figure [4] we show the lowest energy 
levels for the fractional harmonic oscillator for a < 1. This 
spectrum is limited to a finite number of levels, which may 
now be understood in comparison to the plot of zeros of 
the solutions of the free fractional Schrodinger equation 
presented in the lower row of figure [4] 

It can be deduced, that the onset of real eigenvalues 
for a given a of the fractional harmonic oscillator and the 
occurrence of a zero in the potential free solution agree 
qualitatively. 

Since the free fractional Schrodinger equation is for- 
mally equivalent to the classical fractional harmonic os- 
cillator differential equation we may conclude, that in the 
special case of the harmonic oscillator the transition from 
fractional classical mechanics to fractional quantum me- 
chanics is smooth. 

We have thus demonstrated, that the fractional quan- 
tum harmonic oscillator extends the symmetry of the stan- 
dard harmonic oscillator and allows for a generalized de- 
scription of rotational and vibrational degrees of freedom. 

Since we want to emphasize this unifying point of view, 
an appropriate area of application is molecular spectros- 
copy, where an interplay between vibrations and rotations 
is observed. 

Therefore we will investigate the properties of infrared 
spectra of diatomic molecules in the next section using the 
fractional quantum harmonic oscillator. 
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Fig. 4. In the upper row energy levels in the range 0.5 < a < 1 
for the Caputo (left column) and Riemann (right column) frac- 
tional quantum harmonic oscillator are plotted. Blue lines indi- 
cate positive parity, dashed red lines indicate odd parity of the 
corresponding solution. For decreasing a, the number of real 
eigenvalues is reduced and are replaced by a conjugated pair 
of complex eigen values, dotted lines show the real part Re(E) 
and dashed lines the complex Re(E)±Im(E). Green lines show 
the WKB-approximation E WKB . The lower row shows the zeros 
of eigen functions for the potential free fractional Schrodinger 
equation. 



tion are able to perform vibrations around this equilibrium 
position. 

For higher excitations, anharmonic contributions and 
the influence of centrifugal stretching are taken into ac- 
count via a series expansion, which we present here up to 
third order: 



6 Application: The infrared-spectrum of H 3 ^CI E rot _ vih (y,J) - E rot + £ vib - (34) 

£,J(J + l) + A,J 2 (J + l) 2 + 

Traditionally diatomic molecules are treated in lowest or- 111 
der as rigid rotors with a fixed bond length, which in addi- %uj e (v + -) + hoj e x e (v + -) 2 + huj e y e (v + -) 3 
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where the constants B v and D v add a rot-vib coupling 
via: 



0.8 



B v =B e -a e (v+-) 
D u = D e -f3 e (is+^) 



(35) 
(36) 



which results in a seven parameter energy formula to de- 
termine the energy levels in the rot-vib model. 

Now we introduce the fractional quantum harmonic os- 
cillator model, which describes rotational and vibrational 
spectra from a unifying point of view. As we pointed out 
in the previous section, depending on the value of the frac- 
tional derivative parameter a we may generate both types 
of spectra: 



E(a « 2, J) 
E(a w 



(37) 
(38) 



Thus we propose the fractional analogue of the standard 
rot-vib-model: 



B e E(aj « 2, J) + ftCo e E(a v « 1, v) + 
a e £(aj w 2, J)£K » 1^) + £o (39) 

with six parameters aj, B ei a ui Tiuu e , a e , go. A least square 
fit of the fractional rot-vib model energies with the stan- 
dard rot-vib model energies yields the parameter sets listed 
in table H 

In order to emphasize the property of the fractional 
rot-vib model to realize rotational type spectra, we want 
to explore the lowest energy photon absorption processes 
in the infrared region. In that case from A J = ±1, we 
obtain two branches (R and P) with the properties: 



AE R (J)=E(v=l,J+l) 
AEp(J) = E(y = 1, J — 1) 



E[y = 0, J) J = 0, ...(40) 
E(y = 0, J) J = 1,...(41) 



These transition energies may be directly compared to the 
experimental spectrum. The corresponding intensities Ij 
for a given transition are determined in the standard case 
by the Boltzmann-distribution 



I j = (2J+l)e~ 



t/kT 



(2J+l)e-^ J ( J+1 )/ feT (42) 



where k is the Boltzmann constant. The factor (2 J + 1) 
is a result of the multiplicity of a given rotational state. 
In the case of the fractional rot-vib model this implies the 
following generalization: 



7"fr actio 

3 



(2(aj-l)J+l)e- 



t/kT 



(43) 



where the factor (2 (a — 1) J + 1) now may be interpreted 
physically as a smooth transition of multiplicities from 1 
to 2 J + 1 in the range 1 < a < 2, which means the tran- 
sition from U(l) to SO(3). This may be interpreted as an 
indication for the connection of fractal space dimensions 
and fractional calculus. 
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Fig. 5. Fit results for the infrared spectrum of H 35 C1 in the 
wavelength range from 2600-3100 [cm -1 ]. Black lines are the 
experimental values from [16] . From top to bottom the result 
for the energy levels of the fractional rot-vib-model ( 39 ) using 
the WKB-approximation, the Riemann and the Caputo frac- 
tional derivative definition. 



With intensities (43) and assuming an exponential dis- 
tribution for every transition we are able to fit the exper- 
imental spectrum [16] . The results are presented in figure 

El 

The overall error is less than 4%. As a remarkable fact, 
the experimental spectrum is reproduced in the case of 
the Riemann fractional harmonic oscillator with aj > 2, 
which means, the corresponding eigen functions are only 
box-normalizable. Consequently, only the Caputo fractio- 
nal rot-vib model describes the experiment within the al- 
lowed parameter range. 

It is remarkable, that in the case of the fractional rigid 
rotor [17 , which is an alternative description of rot-vib 
systems, just the opposite conclusion was drawn: the Ca- 
puto derivative based version was discarded and the Rie- 
mann rotor model was successful. 
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7 Conclusion and outlook 

It has been demonstrated, that the fractional quantum 
harmonic oscillator generates an energy spectrum, which, 
besides vibrational degrees of freedom also shows rota- 
tional type spectra. The fractional parameter a allows for 
a smooth transition between these two extreme cases. 

Therefore vibrations and rotations may be treated equi- 
valently in a unified generalized rot-vib model, which may 
be successfully used to reproduce the infrared spectra of 
diatomic molecules, which has been demonstrated for the 
case of hydrogen chloride. 

It has also been shown, that the WKB- approximation 
is a useful tool to approximate higher energy levels, but 
in the low energy region the exact solutions using the Rie- 
mann and Caputo derivative definitions differ significantly 
and cover a much broader area of useful applications. 

In this paper we presented first results for the frac- 
tional quantum harmonic oscillator. The results encour- 
age further studies in this field, especially the knowledge 
of eigen values and eigen functions for higher energy lev- 
els [18] will be useful to describe highly excited rotational 
molecular states, an area, which has been made accessible 
recently, see e.g. [19] . 
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